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SECTION  I 


INTRODUCTION 


When  performing  large  finite  element  analysis  it  is  often  desirable 
to  use  a postprocessor  to  identify  those  elements  within  the  structure 
which  are  most  critically  loaded.  The  stress  recovery  package  described 
in  this  report  identifies  those  elements  within  the  structure  having 
safety  margins  less  than  a user  specified  value.  This  package  can  be 
used  as  a part  of  the  initial  NASTRAN  analysis  run  or  interactively  at 
some  later  time. 

This  package  is  composed  of  three  programs.  Program  STRESSA  extracts 
the  stress  data  from  the  NASTRAN  output.  Program  YIELD  uses  the  NASTRAN 
data  deck  to  generate  failure  stress  data  for  membrane  and  plate  elements. 
Program  STRESSB  uses  the  data  generated  by  STRESSA  and  YIELD  to  calculate 
safety  margins  for  the  elements  and  identify  those  elements  with  safety 
margins  less  than  a user  specified  value.  These  programs  are  currently 
operational  on  the  CDC  6000  series  computers.  The  programs  are  coded  in 
FORTRAN  Extended- however,  some  modifications  may  be  required  for  use  with 
other  computer  systems. 


SECTION  II 


The  NASTRAN  Run 

Several  changes  must  be  made  when  running  NASTRAN  with  this  package. 

The  TABPT  module  in  NASTRAN  allows  data  blocks  to  be  stored  in  a compressed 
format  on  the  output  file.  In  order  to  use  this  feature  a call  to  TABPT 
must  be  placed  in  the  DMAP  sequence.  This  call  should  be  after  the  SAVE 
CARDNO  $ command  following  the  call  to  SDR2  (the  stress  recovery  package). 

This  can  be  accomplished  by  placing  an  Alter  Package  in  the  Executive 
Control  Deck  such  as  the  following  (for  Rigid  Format  1,0): 

ALTER  121 

TABPT  0ES1 , ,, , //  $ 

ENDALTER 

NASTRAN  only  calculates  element  stresses  when  stress  output  is  requested 
in  the  Case  Control  Deck.  This  can  be  done  either  with  a STRESS=ALL  or 
STRESS(PUNCH)=ALL  command.  These  commands  allow  the  stresses  to  be  calculated 
with  or  without  printed  stress  output  respectively. 

The  following  NASTRAN  load  and  execute  command  should  be  used:  NASTRAN 
(0ATA1.DATA2, DUMMY) ATTACH.  The  local  file  DATA!  contains  the  NASTRAN  bulk 
data  deck  to  be  executed.  DATA2  will  contain  the  entire  NASTRAN  output  file 
from  which  the  compressed  stress  table  printed  by  TABPT  will  be  removed  in 
STRESSA.  DUMMY  contains  the  stress  data  requested  with  the  STRESS  (PUNCH)=ALL 
card  and  can  be  discarded.  The  file  DATA1  must  be  at  the  beginning  of  infor- 
mation when  NASTRAN  is  executed.  All  the  other  data  files  used  by  the  pro- 
grams in  this  package  except  the  input  file  to  STRESSB  are  rewound  as  required 
internally. 


SECTION  III 
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STRESSA 


The  program  STRESSA  reads  the  NASTRAN  output  file,  removes  the  stress 
table  generated  with  the  TABPT  module,  and  places  the  remainder  on  the 
system  output  file.  A data  file  (TAPE3)  is  generated  containing  the 
stress  information  in  the  format  used  by  STRESSB.  If  the  user  wishes  to 
use  the  stress  recovery  package  interactively  rather  than  as  a part  of  the 
NASTRAN  run  the  program  STRESSA  should  be  executed  after  the  NASTRAN 
execution  and  TAPE3  stored  for  later  use. 
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SECTION  IV 


YIELD 

NASTRAN  does  not  calculate  safety  margins  for  membrane  and  plate 
elements.  Program  STRESSB  (SUBROUTINE  FAILURE)  generates  these  safety 
margins  using  failure  stress  data. 

Program  YIELD  uses  the  NASTRAN  data  deck  to  generate  failure  stress 
data  for  the  membrane  and  plate  elements  based  on  values  given  on  MAT1 
and  MAT2  cards.  These  values  are  then  used  by  program  STRESSB  along  with 
the  stresses  from  STRESSA  to  determine  the  safety  margin  which  is  not 
calculated  bv  NASTRAN  for  these  elements.  YIELD  currently  handles  the 
following  element  types 


CQDMEM 

CTRIA1 

CQUADl 

CTRMEM 

CTRIAZ 

CQUAD2 

CTRBSC 

CTRPLT 

CQDPLT 

Some  requirements  are  placed  on  the  NASTRAN  data  file  for  use  by 
program  YIELD.  Card  identifiers  (col.  1-8)  must  be  left  justified  for 
element  cards,  property  cards,  and  MAT1  and  MAT2  cards.  Exponents  used 
in  yield  stress  values  on  MAT1  and  MAT2  cards  must  be  right  justified. 
Continuation  cards  for  MAT1  and  MAT2  cards  must  immediately  follow  their 
parents. 


SECTION  V 
STRESSB 

Program  STRESSB  lists  all  elements  whose  safety  margins  are  below  a 
user  specified  value.  Output  for  these  elements  consists  of  the  element 
number,  element  type  and  the  stresses  in  that  element.  All  stresses 
normally  output  by  NASTRAN  are  output  for  these  elements.  STRESSB  will 
process  the  following  element  types 


CROD 

CONROD 

CTRIA1 

CQDMEM 

CBEAM 

CELAS1 

CTRBSC 

CTRIA2 

CTUBE 

CELAS2 

CTRPLT 

CQUAD2 

CSHEAR 

CELAS3 

CTRMEM 

CQUAD1 

CTWIST 

CELAS4 

CQDPLT 

CBAR 

Where  possible  STRESSB  uses  safety  margins  calculated  by  NASTRAN.  The 
safety  margins  of  plate  and  membrane  elements  are  determined  in  subroutine 
FAILURE.  This  subroutine  checks  the  major  and  minor  principal  stresses 
against  the  tensile  or  compressive  failure  stresses  according  to  their 
sign.  FAILURE  also  compares  the  maximum  shear  stress  with  the  shear  failure 
stress.  If  any  of  these  are  below  the  user  specified  value  the  element  is 
listed.  For  bending  elements  this  is  done  for  both  surfaces.  The  user  may 
supply  his  own  subroutine  FAILURE.  Elements  not  included  in  the  above  list 
and  elements  for  which  no  failure  stresses  have  been  supplied  will  be 
ignored  by  STRESSB. 

STRESSB  examines  the  element  stresses  by  subcase.  When  the  program 
is  executed  interactively,  safety  margin  and  subcase  number  combinations 


. T.i  -A’/L.' 


are  input  directly  from  the  terminal  when  requested  by  the  program.  When 
STRESSB  is  executed  as  part  of  the  NASTRAN  analysis  run  it  is  necessary  to 
use  an  input  file  with  two  card  images  for  each  subcase  to  be  examined; 
one  card  image  contains  the  real  value  of  the  safety  margin  and  the  integer 
value  of  the  subcase  number  separated  by  a comna,  the  other  card  image  con- 
tains the  word  yes  or  no  indicating  whether  there  is  or  is  not  a subsequent 
subcase  to  be  examined. 

Some  users  may  wish  to  input  the  failure  stress  tables  directly  rather 
than  use  program  YIELD.  This  can  be  done  by  providing  TAPE4  and  TAPE5. 

TAPE4  must  contain  one  card  image  for  each  membrane  or  plate  element  to  be 
examined  containing  the  Element  Number  and  the  Failure  Set  Identification 
Number  of  the  appropriate  Failure  Stress  Set  (215).  TAPE5  must  contain  one 
card  image  for  each  set  of  unique  Failure  Stresses  containing  the  Failure 
Set  Identification  Number  and  the  Failure  Stresses  in  Tension,  Compression, 
and  Shear  (I5.3F13.0). 
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APPENDIX  A 
PROGRAM  STRESSA 


Purpose:  To  rtcover  the  compressed  stress  data  from  a NASTRAN  output 
prior  to  printing. 

External  References:  None 


•DECK 


C 

C 

C 

c 

c 

c 

c 
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10 

11 


101 
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18 

13 

14 

150 

15 

16 
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20 

50 

70 

71 


STRESSA 

PROGRAM  STRESSA< INPUT, OUT  PUT , T APE1  = INPU T, TAPE2 ^OUTPUT , TAPE3 
,T APE6,TAPE7) 

DIMENSION  A(14),C(10>,D(10>,IV(50> 

REWINO  1 
111=0 


THIS  SUBROUTINE  STRIPS  THE 
NASTRAN  OUTPUT  FILE (INPUT) , 

PLACES  THE  STRESS  DATA  ON  A 
DATA  FI LE ( T APE3) , AND  THE 
REMAINDER  OF  THE  NASTRAN  OUTPUT 
ON  THE  SYSTEM  OUTPUT  FILE ( OUTPUT) , 

DIMENSIONS  IN  THIS  PROGRAM 
ARE  NOT  PROBLEM  DEPENDENT 


RE  AD ( 1 » 2 ) A 
FORMAKAl,  13A1C) 

IF (EOF ( 1) ) 1000,3 

IF (A (2) »EQ,1 OH RECORD  NO,)GO  TO  10 
WRITE(2,2>  A 
GO  TO  1 
111  = 1 
IOUT=3 
READ (1,2) A 

IF(A(2) .NE.10HRECORD  NO.)GO  TO  11 
IF(III.EQ.2)GO  TO  150 
READ (1,102) ISK 
FORMAT (Al) 

I F ( E OF ( 1 ) ) 100,103 
IF(ISK.NE.1H  ) GO  TO  101 
BACKSPACE  1 

READ (1,1 2) ITYPE,ISUB,NWDS 
FOf=MAT(2  7X,2I13,65X,I13> 

I F ( I SUB, EQ • 0 ) GO  TO  101 
WRITE (IOUT , 18) ITYPE, ISU8  »NWDS 
FORMAT (10HNEW  RECORD, 3110) 


READ (1,14)3 
FORMAT ( 1 X, A10 ) 

IF(B.NE. 10HRECORD  NO.)GO  TO  13 
IF (III. EQ. 2) WRITE (IOUT, 16) 

RE  AD (1,16) ODD,  (C(I> ,0(1) ,1  = 1,10) 
FORMAT(A1,10(A10,A3) ) 

IF(OOO.EQ.lHl)GO  TO  70 

IF  (C(l) .EQ.IOHRECORO  NO.)GO  TO  50 

OO  20  1=1,10 

IF <C(I) . EQ. 10H  . ANO.O(I) .CQ. 

WRITE (IOUT ,17) C(I),D(I) 

FORMAT (10X ,A10 ,A3) 

CONTINUE 
GO  TO  15 
BACKSPACE  1 
GO  TO  11 
RE  AO (1,71) 

F ORM AT ( 1 OX ) 

REAO (1,71) 


3 H ) GO  TO  15 
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APPENDIX  B 


PROGRAM  YIELD 

Purpose:  This  program  generates  failure  stress  data  for  NASTRAN  membrane 
and  plate  values  from  data  available  on  MAT1  and  MAT2  cards. 

Arrays:  I PROP  (NPROP)  where  NPROP  is  the  largest  membrane  or  plate  property 
number. 

External  References:  RIJUST 
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♦DECK  YIELD 

PROGRAM  YIELD  (INPUT  , OUTPU T , T APE  1= IN» UT , TAPE  2 = 0 UTPUT , TAPE4 , TAPE5, 
IT  APE6) 

DIMENSION  A (12) , ELEM ( 9 ) , IPRO P ( 20 0 0 > 

INTEGER  RIJUST 


IPP.OP  SHOULD  BE  DIMENSIONED  AS  FOLLOWS* 


IPROP ( NPROP)  WHERE  NPROP  = THE  LARGEST  MEMBRANE 
OR  PLATE  PROPERTY  NUMBER 

ELEM(1)=7HQ0MEM 
ELEM(2)=7HTRMEM 
ELEM(3)=7HTRIA1 
ELEM(4)=7HTRBSC 
ELEM(5)=7HTRPLT 
ELEM(6)=7HQ0PLT 
ELEM(7)=7HTRIA2 
ELEM(8)=7HQUA02 
ELEM(9)=7HQUA01 
DO  10  JJ= 1,9 
REWIND  1 

1 REAO(l,2)A 

2 F0RMAT(A1,A7,8A8, A1,A7) 

IF (A(2).EQ.7HN0DATA  )GO  TO  5 
IF (A(2) .NE.ELEM(JJ) ) GO  TO  1 
IF(A(1).EQ.1HP)G0  TO  4 
IELEM=RI JUST (A(3) ) 

IPR01=RI JUST(A(4) ) 

WRITE(6,3> IELEM, IPROl 

3 FORMAT  (215) 

GO  TO  1 

4 I PROP (RIJUST ( A ( 3) )>=RIJUST (A (4)) 

GO  TO  1 

5 REWIND  6 

6 READ(6, 3) IELEM, IPROl 
IF (IELEI.EQ.O)GO  TO  8 
IF  (EOF  (6)  ) 8, 7 

7 WRITE(4,3) IELEM, IPROP ( IPROl ) 

GO  TO  6 

8 REWIND  6 
WRITE  (6,21) 

21  FORMAT ( 10  X ) 

REWINO  6 
00  9 J=l, 2000 

9 IPROP(J)=0 

10  CONTINUE 
REWIND  1 

11  RE AO ( 1, 2 ) A 

IF (A (2) • EQ.7HNDDATA  )GO  TO  20 
IF(A(2).EQ.7HAT1  )GO  TO  12 
IF (A (2 ) • EQ • 7HAT2  ) GO  TO  18 
GO  TO  11 


12  J=  RI JUST ( A ( 3 ) ) 
READ (1,13) TEST 


IF(TEST.NE.A(12) ) GO  TO  16 
BACKSPACE  1 
P E AD  ( 1*  1 4)  X , Y , Z 

14  FORMAT ( 8 X , 3 F 8 • 0) 

17  WRITE(5,15)J,X,Y,7 

15  F0RMAT(I5, 3F13.2) 

GO  TO  11 

16  X = G . 0 
Y=C.O 

z=c.o 

GO  TO  17 

lfi  J=  F I JUST ( A ( 3 ) ) 

READ (1,13) TEST 

IF (TEST  »NE» A ( 12) ) GO  TO  16 

BACKSPACE  1 

RE  AD ( 1, 19)  X , Y » Z 

19  F0RMAT(48X,3F8.C)  i' 

GO  TO  17 

20  STOP 
ENC 


FUNCTION  RIJUST 


Entry  Point:  RIJUST 


Purpose:  To  convert  an  alpha  variable  to  integer  format  with  right 
justification. 


Calling  Sequence:  I1=RIJUST( IWORD) 

IWORD  = An  alpha  format  of  an  integer  variable 
II  - The  right  justified  integer 


INTEGER  FUNCTION  RI JUST< IWORO) 
M A SK  = 553 
00  10  1=1,10 
ITEST=77B#AND#I WORD 
IF(MASK.NE.ITEST)GO  TO  2 C 
10  IW  OR  D= SHIFT  (IWCRD»54) 

20  OECOOEdO  * 30  , IWORD)RIJUST 

30  FORHATdlOl 

RETURN 
END 
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PROGRAM  STRESSB 

Purpose:  This  program  examines  the  stress  data  from  a NASTRAN  execution 
and  prints  out  those  elements  for  which  the  safety  margin  is  below  a user 
specified  value. 

Blank  Common: 

I SAFT ( NSAFT ) , SAFMAR , SAFTVAL ( NVAL 3 ) , X ( 20 ) 

NSAFT  = The  largest  membrane  or  plate  element  number. 

ISAFT(NN)  = The  property  number  of  plate  or  membrane  element  NN. 

NVAL  = The  largest  property  number  of  a plate  or  membrane  element. 
SAFTVAL(NN.J)  = The  J'th  failure  stress  for  elements  with  property 
number  NN  (tension,  compression,  and  shear). 

SAFMAR  = The  user  input  critical  safety  margin. 

X(J)  - The  stresses  in  the  current  element. 

External  References:  CROD,CSHEAR,CMEM,ELAS,CBEAM 
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♦OECK  STRESSB 

PROGRAM  STRESSB ( I NPUT, OUTPUT, TAPE  1 = 1 NPUT. TAPE  2=0 UTF,. , , 
1TAPE3,TAPE4,TAPE5> 

COMMON  ISAFT<  100001  , S ARM AR, S AFT V AL ( 1 C 0 0 , 3) , X(20) 

THE  VARIABLES  SHOULD  BE  DIMENSIONED  AS  FOLLOWS 

ISAC  T < NSAFT ) 

SAFTVAL(N\/AL,3> 

NSAFT  = THE  LARGEST  MEMBRANE  OR  PLATE  ELEMENT  NUMBER 
NVAL=THE  LARGEST  MEMBRANE  OR  PLATE  PROPERTY  NUMBER 


4 

5 

1 


REWIND  1 
REWIND  4 

R E AD (4,2) I, ISAFT(I) 

FORMAT (215) 

IF  (EOF (4)  ) 33,1 
REWIND  5 

READ(5,35> I, (SAFTVAL(I,J) »J=1,3) 

FORMAT (15, 3E13.6) 

IF  (£  OF ( 5 ) ) 31,34 
WRITE(2,4) 

FORMAT (”  WHAT  IS  THE  CRITICAL  SAFETY 
READU,*)  SAFMAR,INUMB 
REWIND  3 
READ (3,6) A 
FORMAT (A10 , 3110) 

IF (A.NE. 1QHNEW  RECORD ) GO  TO  5 
BACKSPACE  3 

REA0(3,6) A,ITYPE,ISUB,NWOS 
IF(ITYPE.EO.O)GO  TO  100 
IF(ISUB.NE.INUMB)GO  TO  5 
IF(ITYPE.EQ.l) CALL  CROC 
IF(ITYPE.EQ.2)CALL  CBEAM(ITYPE) 

IF  (ITYPE.EQ.3) CALL  CROD 
IF  (ITYPE.EQ.4) CALL  CSHEAR 
IF(ITYPE.EQ.5) CALL  CSHEAR 
IF(ITYPE.EQ.6)CALL  CPLATE 
IF(ITYPE.EQ.T) CALL  CPLATE 
IF  (ITYPE. EO. 8) CALL  CPLATE 
IF  (ITYPE. EQ. 9) CALL  CME M 
IF  (ITYPE .EQ.10) CALL  CROD 
IF  (ITYPE.EQ.ll)CALL  ELAS 
IF  (ITYPE,EQ,12) CALL  ELAS 
IF(ITYPE.EQ.13)CALL  ELAS 
IF(ITYPE.EQ.14)GO  TO  5 
IFdTYPE.EQ.15)  CALL  CPLATE 
IF  (ITYPE .EQ.16) CALL  CMEM 
IF  (ITYPE. E0.17)CALL  CPLATE 
IF(ITYPE.E0.16)CALL  CPLATE 
IF(ITYPE.EQ.19)CALL  CPLATE 
IF  (ITYPE. GT .19.ANn.ITYPE.LT. 34)G0  TO  5 
IF  (ITYPE. EQ.34)CALL  CBEAM(ITYPE) 
IF(ITYPE.GT.34)G0  TO  «? 

GO  TO  6 
WRITE  < 2, 10  1) 


16 


MARGIN  AND  SUBCASE  NUMBER?") 


# 

v#' 

V 


130 


I 


SUBROUTINE  CROD 


Entry  Point:  CROD 


Purpose:  To  determine  if  a rod's  element  stresses  are  below  the  critical 
value. 


Calling  Sequence: 
CALL  CROD 


Blank  Common:  See  main  program 
X(J)  not  used 


v 


SUBROUTINE  CROO 

COMMON  IS AFT ( 10  000) ,SAFMAR,S AFT VAL (1000 ,3) 

WRIT  E (2 » 3 ) 

FORMAT  ( 1H1 , " ROD  ELEMENTS  **/“  STRESS  SAFETY  MARGIN") 

READ(3,1) A.IELEM 

FORMAT(A10,I12) 

IF(A,EQ.10HNEW  RECORO)GO  TO  10 
READ  13.2) AXIAL 
FORMAT(10X,F13.0) 

RE  AO  (3,2) ASM 
IF  (ASM.EQ.l.O)GO  TO  20 
IF  (ASM. GT.SAFMAR)GO  TO  20 
WR ITE  (2 , U ) IELEM, AXIAL, ASM 
FOKMAT(1X,HO,E13.6,F13.  8,"  AXIAL") 

READ<3,2) AXIAL 

READ(3,2) ASM 

IF  (ASM.EQ. 1. 3)GO  TO  6 

IF (ASM.GT .SAFMAR)GO  TO  6 

WRITE (2, 7) IELEM, AXIAL, ASM 

FORMAT (IX, 1 10, El 3. 6, FI 3. 8 ,"  TORSIONAL") 

GO  TO  6 
BACKSPACE  3 
RETURN 
END 


OK' 
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SUBROUTINE  CSHEAR 


Entry  Point:  CSHEAR 

Purpose:  To  determine  if  the  stress  in  a shear  panel  is  below  the 
critical  level. 

Calling  Sequence 
CALL  CSHEAR 

Blank  Common:  See  main  program 
X(l)  = Maximum  Shear  Stress 
X(2)  = Average  Shear  Stress 
X(3)  = Safety  Margin 
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SUBROUTINE  CSHEAR 

COMMON  ISAFT(1C00G),SAFMAR,SAFTVAL(1CCG,3),X(3> 
WRITE<2, 3) 

3 FORMAT  (1H1  ,**  SHEAR  PANELS’*/ 

1"  ELEMENT  MAX  SHEAR  A VE  SHEAR  SAFETY  MARGIN”) 

1 REAOI.3,2)  A.IELEM 

2 FORMAT ( A 1 0 , II 2) 

IF(A.EQ.  1CHNEW  RECOROGO  TO  10 

00  ft  1=1,3 

READ(3,4)X(I) 

A FORMAT (10X,F13.0) 

CONTINUE 

IF(X(3).EO.1.0)GO  TO  1 
IF(X(3) ,GT.SAFMAF)GO  TO  1 
WRIT  E ( 2 , 7 ) IELEM,  X 
FORMAT  <16* 3E13.6) 

GO  TO  1 
BACKSPACE  3 
RE  TURN 


I 

SUBROUTINE  CMEM 

Entry  Point:  CMEM 

Purpose:  To  determine  if  the  stress  field  in  a membrane  element  is 
critical. 

External  References:  FAILURE 

Calling  Sequence: 

CALL  CMEM 

Blank  Common:  See  Main  Program 
X(l)  = not  used 
X(2)  = Norma 1-X 
X(3)  = Normal-Y 
X(4)  = Shear-XY 
X(5)  = 0-Shear  Angle 
X(6)  = Major  Principal 
X(7)  = Minor  Principal 
X(8)  = Maximum  Shear 


SUBROUTINE  CMEM 

COMMON  ISAFT(  10  0 0 0) , S AFM AR , S AFT V AL (1000 ,3) , X C 8) 
WRITE  (2*  3) 

3 FORMAT ( 1H1 , " MEMBRANE  ELEMENTS"/"  ELEMENT  NUMBER"/ 

117X,"N0RMAL-X",6X,"N0RMAL-Y",6X,“SHEAR-XY",7X, 

2” A NCLE“, 7X, "MA JPRIN", 6X , "M IN PRIN",2X,“ MAX-SHEAR") 

1 READ<3,2> A.IELEM 

2 FORMAT (A10 ,112) 

IF ( A . EO. 10 HNE W RECORD) GO  TO  10 
OO  6 1=2,8 
READ ( 3,4 ) X ( I) 

A FORMAT (10X,F13,0) 

6 CONTINUE 

IT  = 1 

CALL  FAILURE(IELEM, IT) 

GO  TO  1 

10  BACKSPACE  3 

RETURN 
END 
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SUBROUTINE  FAILURE 


Entry  Point:  FAILURE 

Purpose:  To  calculate  safety  margins  for  plate  and  membrane  elements. 

Calling  Sequence: 

CALL  FAILURE  (IELEM.IT) 

IELEM  = the  element  number 
IT  = 1 for  membrane  elements 
= 2 for  plate  elements 

Blank  Common:  Sqe  main  program 
X(J)  see  CMEM  or  CPLATE 


oooooooooaoooooooooooooooonooooonoooo 
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SUBROUTINE  FAILURE (IELEM.IT) 

COMMON  ISAFT  ( 1C0C0)  , SAFM  AR , S AFT  V AL  < 1 00  G , 3)  , X < 16) 


THIS  SUBROUTINE  OEFINES  THE  FAILURE  CRITERIA 
FOR  PLATE  AND  MEMBRANE  ELEMENTS 


THE  USER  MAY  REPLACE  THE  CODE  BETWEEN  HERE  ANO  STATEMENT 
20  WITH  HIS  OWN  ALGORITHM  TO  EXAMINE  THE  STRESS 
FIELD  ANO  OETERMINE  WHETHER  TO  PRINT  OUT  THE  STRESSES 
FOR  THE  ELEMENT (GO  TO  20)  OR  NOTIGO  TO  30).  THE  FOL- 
LOWING VARIABLES  ARE  AVAILABLE  FOR  USE 


1 


7 

3 


IELEM=ELEMENT  NUMBER 
IT=1  FOR  MEMBRANE  ELEMENTS 
IT=2  FOR  PLATE  ELEMENTS 

ISAFT <IELEM)=THE  FAILURE  SET  10  FOR  ELEMENT  IELEM 
SAFMAR  =THE  USER  SPEC  IFIED  C RITIC  AL  SAFETY  MARGIN 


SAFTVALdSAFT(IELEM)  ,1)  = TENSI LE  FAILURE  STRESS 

2=C0MPRESSIVE  FAILURE  STRESS 
3=  SHEAR  FAILURE  STRESS 


MEMBRANES 

PLATES 

X(l) 

0.0 

Zl  FIBRE  DISTANCE  1 

2 

NORMAL-X 

NORMAL-X  AT  Z1 

3 

NORMAL-Y 

NORMAL-Y  AT  Z1 

4 

SHEAR-XY 

SHEAR-XY  AT  Zi 

5 

0-SHEAR  ANGLE 

0-SHEAR  ANGLE  AT  Zl 

6 

MAJ.  PRIN. 

MAJ. PRIN.  AT  Zl 

7 

MIN.  PRIN. 

MIN. PRIN.  AT  Zl 

8 

MAX.  SHEAR 

MAX. SHEAR  AT  Zl 

9 

NOT  USED 

72  FIBRE  DISTANCE  2 

10 

NOT  USED 

NORMAL-X  AT  Z2 

11 

NOT  USED 

NORMAL-Y  AT  Z2 

12 

NOT  USED 

SHEAR-XY  AT  Z2 

13 

NOT  USEO 

0-SHEAR  ANGLE  AT  Z 2 

14 

NOT  USEO 

MAJ. PRIN.  AT  Z2 

15 

NOT  USED 

MIN. PRIN.  AT  Z2 

16 

NOT  USED 

MAX  . SHEAR  AT  Z2 

SAFT  C 

IELEM)  , EQ. 0 > RETURN 

SAFETYT=SAFTVAL< ISAFT (IELEM)  ,1) 
SAFETYC=-SAFTVAL( ISAFT (IELEM) , 2) 
SAFETYS=SAFTVAL(ISAFT(I£LEM)  ,3) 

N = 1 


NN  = N«-5 


NNN=N«-6 


□ 0 3 1=  NN » NNN 
IF(X(I).EQ.O)GO  TO  3 
IF  (X (I) ,GT .0.0)GO  TO  7 
SAFE=SAFETYC/X(I) 

IF  (SAFE.LE.SAFMAR)GO  TO  20 
GO  TO  3 

SAFE=SAFETYT/X(1) 

IF  (SAFE.LE.SAFMAR)GO  TO  20 
CONTINUE 

IF(X(I).E0.3)RETURN 
NN=N  *7 
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SAFE=SAFETYS/A3S(X(I)  ) 
IF(SAFE.LE.SAFMAR)GO  TO  20 
IF  (IT.EQ  .1 ) RETURN 
IF  (N.EQ.9) RETURN 
N=9 

GO  TO  1 

WF IT  E ( 2*  22 ) I ELEM 
FORMAT(1Hj*I10) 

WR  IT  E ( 2 » 2 1 ) <XCI)  ,1=1,8) 
FORMAT (1X,6E13.6) 

IF  (IT .E3.1) RETURN 

WRITE (2, 21)  <X(  I)  , 1=9, 16) 

RETURN 

ENO 
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SUBROUTINE  ELAS 


Entry  Point:  ELAS 

Purpose:  To  determine  if  the  stress  in  an  elastic  spring  element  has 
reached  the  critical  level. 

Calling  Sequence: 

CALL  ELAS 

Blank  Common:  See  Main  Program 
X(J)  - not  used 


SUBROUTINE  ELAS 

COMMON  ISAFT ( 10000) , S AFM AR, S AFT VAL ( 1 0 0 0 ,3) 

WRITE (2*3) 

FORMAT ( 1H 1 1 **  ELASTIC  SPRINGS"/**  ELEMENT  STRESS**) 

REAO(3»2) A.IELEM 

FOFMAT(A10,I12) 

IFIA.EQ. 10HNEW  RECORD) GO  TO  1G 

RE  AD  ( 3*  4 ) X 

FORMAT(10X,F13.C) 

IF  (X .GT. 0 • 0 ) GO  TO  5 

IF  (ISAFT (IELEM) . EQ.O .0 )GO  TO  1 

S A FETY=SAFT VAL (ISAFT (IELEM)  ,2) 

SAFE=SAFETY/X 
IF(SAFE.GT.SAFMAR)GO  TO  1 
WRITE(2,6> IELEM, X 
FORMAT (IX  tl/t£l3>6) 

GO  TO  1 

IFdSAFT  (IELEM)  .EQ«G«0)GO  TO  1 
S AFETY  = S AFT VAL (ISAFT (IEuEM»  , i> 

SAFE= (SAFE TY-X) /SAFETY 

IF  (SAFE.LE.SAFMAR)GO  TO  7 

GO  TO  1 

BACKSPACE  3 

RETURN 

END 


# 
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SUBROUTINE  CPLATE 


Entry  Point:  CPLATE 

Purpose:  To  determine  if  the  stress  in  a plate  element  has  reached  the 
critical  level. 

Calling  Sequence: 

CALL  CPLATE 

External  Reference:  FAILURE 


Blank  Common: 

n 

See  Main 

12 

Program 

X(l) 

X ( 9) 

Fibre  Distance 

X(2) 

X(10) 

Norma 1-X 

X ( 3) 

X(11) 

Normal -Y 

X(4) 

X ( 1 2 ) 

Shear-XY 

X(5) 

X ( 1 3 ) 

0-Shear  Angle 

X(6) 

X ( 1 4) 

Major  Principal 

X(7) 

X ( 1 5 ) 

Minor  Principal 

X(8) 

X(  1 6) 

Maximum  Shear 
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SUBROUTINE  CPLATE 

COMMON  ISAFT<  1C  00  0)  , SAFM AR, SAFTVAL ( 1000,3)  , X<  16) 

XRIT E (2  * 3 ) 

3 FORMAT  ( 1H 1 , " PLATE  ELEMENTS"/”  ELEMENT  NUM3ER"/ 

12  X," FIBER  DISTANCE" ,5 X, "NORM AL-X",6X," NORMAL- Y", 6 X,"SH EAR- XY", 
27X,"ANGLE",7X,"MAJPRIN", 6X ,"MINPRIN", 2 X , "MAX -SHEAR") 

1 READ(3»2)A»IELEM 

2 FORMAT(A10, 112) 

IF  (A  .EQ. 10HNEM  RECORD) GO  TO  10 
DO  6 1=1,16 
REA0(3,4)X (I) 

FOrMAT(10X»Fl3.0) 

CONTINUE 
I T=2 

CALL  FAILURE (IELEM, IT) 

GO  TO  1 

11  BACKSPACE  3 

RETURN 

END 


SUBROUTINE  CBEAM 


Entry  Point:  CBEAM 


Purpose:  To  determine  if  a BAR  or  BEAM  element  stresses  have  exceeded 
the  critical  level. 

Calling  Sequence: 

CALL  CBEAM (I TYPE) 

I TYPE  = 2 for  a CBEAM  element 
34  for  a CBAR  element 


Blank  Common: 

See  Main  Program 

X(l) 

SA1 

X ( 9) 

SB1 

X ( 2 ) 

SA2 

x(10) 

SB2 

X ( 3) 

SA3 

X(ll) 

SB3 

X(4) 

SA4 

X ( 1 2) 

SB4 

X ( 5 ) 

Axial 

X ( 1 3 ) 

SB-max 

X(6) 

SA-max 

X ( 1 4 ) 

SB-min 

X(7) 

SA-mi n 

X ( 1 5 ) 

Safety  Margin  in  Compression 

X(8) 

Safety  Margin  in 

Tension 

■i  •£>  *K 
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SUBROUTINE  CBEAM ( IT YPE ) 

COMMON  ISAFT (10000) , SAFMAR.SAFT VAL < 1000 ,3) , X< 15) 

DO  33  1=1,15 
X(I)=G.C 
WRIT  E ( 2, 3) 

FORMAT  (1H1,”  BEAM  ELEMENTS**/"  ELEMENT  NUMBER"/ 
16X,~SAl",10X,"SA2",10X,"SA3",i0X,"SA4",  3X, "AXIAL" , 
28X,"SA-MAX",8X*"SA-MlN"»2X»"SAFETY  MARGIN  IN  TENSION"/ 
36X,"S81”,1QX,"SB2",10X,**SB3",10X,"S94",  8X,"S9-MAX",8X, 
4"SB-MIN**,2X, "SAFETY  MARGIN  IN  COM  PRESS  I ON“> 

READ<3,2> A.IELEM 
FORMAT (A10 , 112) 

IFU.EQ.10HNEW  RECORD)GO  TO  10 
DO  6 1=1,15 

IF(I,EQ,12.AN0.ITYPE.EQ.2)G0  TO  6 
IF(I,EQ,4.AND.ITYPE.E0.2)G0  TO  b 
READ ( 3, 4 ) X ( I > 

FORMAT(10X,F13.0) 

CONTINUE 

IF(X(8).EQ.1.0)GO  TO  15 
IF  < X < 15)  . E 0. 1 . 0 ) GO  TO  20 

IF(X(8) ,GT.SAFMAR.AND.X(15).GT.SAFMAR)G0  TO  1 
WP  ITE (2,7)  IELEM, X 

FORMAT{1HO,I10/1X,8E13.6/1X,7E13.6) 

GO  TO  1 

IF(X(15) ,GT, SAFMAR) GO  TO  1 
IF  ( X ( 15)  . EQ  . 1 • 0)  GO  TO  1 
GO  TO  8 

IF  (X  (8)  ,GT.SAFMAR)GO  TO  1 

GO  TO  8 

BACKSPACE  3 

RETURN 

END 


\ 
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APPENDIX  D 

SAMPLE  CONTROL  CARD  RECORDS 


In  order  to  use  this  package  in  conjunction  with  a NASTRAN  analysis 
run  with  all  the  data  on  the  input  file  the  following  control  cards  are 
required: 

JOBCARD 

ATTACH,  NASTRAN,  **• 

COPYCR,  INPUT,  DATA1 . 

COPYCR,  INPUT,  DATA2. 

NASTRAN  (DATA1 , DATA3,  DUMMY)  ATTACH 
ATTACH,  NASLIB,  **• 

LIBRARY,  NASLIB. 

STRESSA,  DATA3. 

YIELD,  DATA1 . 

STRESSB,  DATA2. 

7/8/9 

-NASTRAN  DATA  DECK  - 
7/8/9 

-STRESSB  DATA  DECK  - 
6/  7 / 8/  9 

Either  DATA1  or  DATA2  could  be  submitted  from  permanent  file/ tape  by 
replacing  the  copy  statements  with  appropriate  ATTACH/REQUEST  statements. 
If  Failure  Stress  Data  is  to  be  input  from  existing  files  the  YIELD,  DATA1 
card  can  be  replaced  with  ATTACH/REQUEST  or  COPY  statements  to  provide 
TAPE4  and  TAPE5. 


m; 


To  generate  a data  file  for  later  interactive  use  in  STRESSB  the 
last  three  statements  of  the  control  record  should  be  replaced  with. 

REQUEST,  TAPE3,  *PF. 

STRESSA. 

CATALOG,  TAPE3,  tape3pfn,  RP=999. 

REQUEST,  TAPE4,  *PF. 

REQUEST,  TAPES,  *PF. 

YIELD,  DATA1 . 

CATALOG,  TAPE4,  tape4pfn,  RP=999. 

CATALOG,  TAPE5,  tapeSpfn,  RP=999. 

If  TAPE4  and  TAPES  are  already  available  for  this  model  the  last  5 statements 
above  can  be  omitted. 
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